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Abstract 

We calculate the vector meson masses in iVf = 2 + 1 Wilson chiral perturbation theory at next-to-leading 
order. Generalizing the framework of heavy vector meson chiral perturbation theory, the quark mass and 
the lattice cutoff dependence of the vector meson masses is derived. Our chiral order counting assumes that 
the lattice cut-off artifacts are of the order of the typical pion momenta, p ~ aAq CD . This counting scheme 
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is consistent with the one in the pseudo scalar meson sector where the 0(a 2 ) terms are included in the 
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■ leading order chiral Lagrangian. 
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I. INTRODUCTION 



This is the second in a series of papers where we compute a variety of mesonic quantities in 2+1 
flavor Wilson Chiral Perturbation Theory (Wilson XPT) . After having computed the pseudo scalar 
mesons masses lj we present here the results for the vector meson masses. The calculation of the 
pseudo scalar decay constants and the axial vector Ward identity quark mass is in progress [jj. 
The main goal of this series of papers is to provide the necessary chiral fit forms for unquenched 
2+1 flavor Lattice QCD simulations with improved Wilson fermions, as they have been currently 
performed by the CP-PACS/JLQCD collaboration Q. 

The quark masses in the CP-PACS/JLQCD simulations are heavier than their physical values. 
The ratio of the pseudo scalar to vector meson mass is in the range mps/my — 0.62 — 0.78. A 
chiral extrapolation in the light up and down quark masses is thus required. In order to perform 
the chiral extrapolation before taking the continuum limit we formulate XPT at non-zero lattice 
spacing, as originally proposed in Refs. Q, Q|. A variety of pseudo scalar quantities has been 
already computed, mainly for 2 flavor Wilson XPT (see Ref. P| and references therein). Although 
the vector meson masses were calculated recently in 2 flavor partially quenched Wilson XPT [7j, 
the result for 2+1 flavors was missing. 

In this paper we follow the heavy vector meson formalism first introduced by Jenkins et. al. 
in Ref. |8jl. The generalization to Lattice QCD at non-zero lattice spacing a is straightforward 

n n 

and mirrors the strategy spelled out in Refs. [J, |a|. However, since the lattice spacing is an 
additional expansion parameter the power counting requires some care. Here we adopt a power 
counting which assumes that lattice cut-off artifacts are of the order of the typical pion momenta, 
aAq CD ~ p, as it was assumed in Ref. Previous results in unquenched 2-flavor simulations [lol ] 
seem to indicate that this is the appropriate power counting for describing the lattice results of the 
CP-PACS/JLQCD collaboration. 

There is a price to pay if one wants to perform the chiral extrapolation before taking the con- 
tinuum limit. The chiral fit forms contain terms proportional to powers of the lattice spacing 
accompanied by additional unknown low-energy constants. These constants are essentially un- 
constrained and serve as additional fit parameters. Obviously, the presence of too many of these 
additional parameters would limit or even spoil the chiral extrapolation. Our one loop results 
for the p and K* meson masses contain seven unknown fit parameters compared to three in the 
corresponding continuum XPT result of Ref. This number is still small enough for the results 
to be useful for the chiral extrapolation of the CP-PACS/JLQCD collaboration data. 
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This paper is organized as follows. In Sec. [Hj we briefly review the heavy meson formalism, as 
it was introduced in Ref. 0. In Sec. 11111 we first summarize the meson and spurion fields which 
are necessary in our calculation. After discussing the power counting we derive the chiral effective 
Lagrangian and compute the vector meson masses to one loop. Concluding remarks are given in 
section II VI 



II. HEAVY VECTOR MESON EFFECTIVE THEORY 



We adopt the heavy vector meson effective theory to derive the quark mass and lattice spac- 
ing dependence of the vector meson masses. 1 This effective theory can deal with vector number 
conserving decay processes like V — > V'X, where V and V are vector mesons, and X represents 
a state with one or more low momentum pseudo scalar mesons. One the other hand, it cannot 
be applied to the process p — > irir, for example, since this decay includes hard pions in the final 
state. 

Some physical quantities have been calculated within the continuum formulation of the heavy 
vector meson effective theory. The first result for the vector meson masses up to 0(p 3 ) can be found 
in Ref. M. The effects of isospin breaking and electromagnetic corrections are included in Ref. 

n n 

121]. The calculation of the vector meson masses to 0(p 4 ) was done in Ref. |l3| . while the results 
for the decay constants are given in Ref. ^| . For results in quenched and partially quenched XPT, 
see Ref. Q] and Ref. [ljj, respectively. 

Let us briefly review the main idea behind the heavy vector meson formalism. The following 
argument can be found in Ref. ^J. The vector meson sector introduces a typical energy scale, 
the vector meson mass my. This mass does not vanish in the chiral limit, and it is of about the 
same size as the chiral symmetry breaking scale, <iirf n ~ 1 GeV. The vector mesons are described 
by heavy matter fields. In order to obtain the Lagrangian of the effective theory one expands the 
relativistic Lagrangian in powers of \jmy. The starting point is the free relativistic vector meson 
Lagrangian, 

Cr = -\v^V, v + \m 2 v V^ (1) 

with Vfn, = dfjVy — dyV^. The relativistic field is decomposed into two parts, a parallel and a 
perpendicular component with respect to the 4-velocity u„ {v 2 = 1) of the vector meson, 

V? = P»vV u +Vn(v ■ V) = V ±li + « M Vf, , (2) 



This formalism is very similar to heavy baryon chiral perturbation theory [Tlj . 



where 



P 



(3) 



is the projector which selects the perpendicular component of the vector field ( P 2 = P, v^P^ u = 0). 



V±n and V\\ can further be written as 

Vu 



\J2my 
1 



[ e -imvvx w ^ + e imvv-x w ^ 



(4) 
(5) 



\/2m v 

where and Wl are effective vector meson fields, and P^vW 11 = and P^W^ = are 
understood. By neglecting terms proportional to exp (±2imyv • x), which oscillate rapidly if one 
takes my to infinity, the Lagrangian Cr in terms of the W fields is given by 



Cr 



1 



2my - 

1 r 
+ 2 



-d^wtd^Wu + (v • aw,f)(u • dWu 



1 



2my . 
iWliy ■ 9)W 



dyWKv ■ dW^) + (v ■ dWhd^Wu 



(6) 



In the infinite mass limit, my — > oo, the W\\ field decouples because of the presence of the mass 
term my W||W||/2. In order to remove the parallel component, one can impose the constraint 
v -V = or v -W = 0. Alternatively, we can remove W\\ by making use of the equation of motion, 



Ww 



—i 



my 



d^ + Oil/ml), 



(7) 



which also shows that W\\ is suppressed by powers of 1/my relative to the perpendicular component 
Wa- Having removed the parallel component the resulting effective Lagrangian simplifies to 



Cr 



-iWUv ■ d)W' 1 + 0(l/m v ). 



(8) 



The first term is the kinetic term for the vector meson in the heavy effective theory. Vector meson 
fields which appear in the following are the effective fields Wn and W^. The partial derivative id^ 
acting on these effective fields produces a small residual momentum r M , defined as 

kf, = myv^ + r M , (9) 

where k^ is the usual four-momentum of the vector meson. We assume here that the residual 
momentum is of the size of the low momentum of the pseudo scalar meson. This assumption holds 
in interaction processes with soft pions. 



III. WILSON CHIRAL PERTURBATION THEORY FOR VECTOR MESONS 



The chiral effective Lagrangian for heavy vector mesons is expanded in powers of the small 
pseudo scalar momenta and masses, and the residual momentum of the vector mesons. The 
symmetries of QCD, in particular chiral symmetry, dictate the form of the terms in the chiral 
Lagrangian. Even though the quark masses explicitly break chiral symmetry, their effect is properly 
taken into account by a spurion analysis, where the quark mass matrix is assumed to transform 
non-trivially under chiral transformations in an intermediate step. 

The same principles apply to Wilson XPT, the low energy effective theory for Lattice QCD with 
Wilson fermions. The main difference is the presence of an additional expansion parameter, the 
lattice spacing a P, ^| . The non-zero lattice spacing is also taken into account by the spurion 
analysis, and the transformation behaviour of the corresponding spurion field is exactly the same 
as for the quark masses. Constructing the terms in the chiral Lagrangian of Wilson XPT for vector 
mesons is therefore straightforward. Apart from the presence of two instead of one spurion field 
there is essentially no difference to Ref. The order counting, however, is not entirely obvious, 
since both the quark masses and the lattice spacing are expansion parameters and their relative 
size is important for a consistent power counting, as we will discuss in section IIIIBI 



A. Matter and spurion fields 

The pseudo scalar meson field is introduced as usual as an SU(3) unitary matrix 

'2ilT 



where 



S = exp 



+ 



f 



(10) 



n = V27r a T a 



ii 



V2 Ve 
K~ 



7T 



K+ 



(11) 



The SU(3) generators T a (a = 1, . . . , 8) are normalized such that 



tv[T a T b ] = -5 a 



(12) 



and / is the leading order pseudo scalar decay constant in the chiral limit. 2 The field £ transforms 
according to S — > LSi? -1 under chiral rotations with L £ SU(3)l and R £ SU(3)r. The square 



Our normalization corresponds to / ~ 132 MeV. 
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root of S, 



cxp 



in 
7 



(13) 



is needed to describe the interaction of the pseudo scalars with the vector mesons. It transforms 
under chiral transformations according to 



U£R-\ 



(14) 



with an SU(3) matrix U. In fact, eq. (|14|1 defines U, which is a function of L, R, and II. For vector 
transformations with L = R one finds U = L = R. 
The vector meson fields are introduced as an octet 



0„ = V2plT a 



9+ K* 

oO ,(8) 



-|- X 



A* 



2d 



(8) 



(15) 



and a singlet = <p° u - These fields are required to satisfy the constraints 



v . S = v ■ O = 0, 



(16) 



in order to describe spin 1 particles with three polarization states. As stated in the previous section, 
this constraint can be enforced by applying the projector P^ u on the vector meson fields. In the 
following we assume that this projector has been applied, i.e. we implicitly assume — Py^S v 



and O m = P^O". 



For the construction of the chiral Lagrangian the quark masses and the lattice spacing are 
treated as spurion fields. For the quark mass matrix we use 



M q = diag(m, m,rh s ) = M I + M 8 T 8 , 

where Mq and Ms are expressed in terms of rh and rh s 

2rh + m s ^ 2{m — m s ) 



(17) 



(18) 



3 ' V3 ' 

We use the tilde in order to highlight that the quark masses rh and rh s denote shifted quark masses, 
which are defined such that the tree level pseudo scalar meson masses become zero if rh = rh s = 0. 
See Ref. and also appendix ^] for details. 

We introduce a spurion field A to include the effect of a non-zero lattice spacing a. This field 
transforms under chiral transformations just as the quark mass field M (M = M q after the spurion 
analysis), i.e. 

-i 



.4 



LAK 



(19) 
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LAB) 
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At 


W± 


uw±W 


wl 


±W± 


v. 




-V T 

A* 


v. 


A, 


UA^W 


A T 


-A, 






-(X>„0„) T 





TABLE I: Transformation properties under the group G = SU{3) L x SU(3) R , (L £ SU(3) L and R e 
SU(3)r), charge conjugation C and parity P. In the parity transformed expressions it is understood that 
the argument is (— x,t). 



Once the chiral Lagrangian is derived, the spurion is set to al where / denotes the unit matrix in 
flavor space. 

For the construction of the chiral Lagrangian it will also be useful to introduce the following 
quantities: 

M± = i(£Mt£±etM£t), (20) 

W± = i(£4t£±£U£t), (21) 

V, = ^d^ + tfdrS), (22) 

A, = (23) 
V u O» = d„0^ + [V v ,0^]. (24) 

The transformation behavior of these quantities as well as of the meson and spurion fields under 
chiral rotations in G = SU(3)l x SU(3)r, charge conjugation C and parity P are summarized in 
table HI 
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B. Power counting 



Having defined the meson and spurion fields it is straightforward to write down the most gen- 
eral chiral Lagrangian which is compatible with chiral symmetry, charge conjugation and parity. 
However, since we have two sources for explicit chiral symmetry breaking, the quark masses and 
the lattice spacing, their relative size matters for an appropriate power counting. Moreover, the 
power counting one wants to employ for the vector meson chiral Lagrangian should be consistent 
with the one adopted for the pure pseudo scalar chiral Lagrangian. 

In Ref. Q] the 0(a 2 ) terms are included in the leading order chiral Lagrangian. To be consistent 
we adopt here the following power counting scheme: 

LO : 0(p),0(a), 

NLO : 0(m),0{a 2 ),0(ap),0(p 2 ), (25) 
NNLO : 0{mp),0(p 3 ),0{am),0{a 3 ),0{ap 2 ),0{a 2 p). 

Here p represents both the residual momentum of the vector meson and the momentum of the 
pseudo scalar meson, and m denotes the quark mass. This counting scheme assumes p ~ a. This 
implies p 2 ~ a 2 , and therefore maintains consistency with the power counting of Ref. Q] for the 
pure pseudo scalar sector. 

We have listed the NNLO contributions in eq. (|25|). because some NNLO terms enter already 
the one-loop calculation of the vector meson masses. To discuss this point, let us define a generic 
power counting scale E of order p ~ a, so that the LO terms or of 0(E), the NLO terms are 
0(E 2 ) and so on. Just from a dimensional analysis one finds that a one-loop integral involving a 
vertex from the LO Lagrangian gives an 0(E 3 ) contribution to the vector meson mass. Therefore, 
some of the coefficients in the NNLO Lagrangian are needed as counterterms for the cancellation 
of the one-loop divergences. The coefficients of the NLO Lagrangian, on the other hand, are not 
needed as counterterms, since there are no one loop corrections with dimension E 2 . This situation 
is different from the pure pseudo scalar sector where the divergences in the one-loop contribution 
from the LO Lagrangian are canceled by the tree level terms coming from the NLO Lagrangian. 
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The structure of the correction to the vector meson masses can be summarized as follows: 



my = n c + [tree level contribution from the LO Lagrangian] 

^ „ ? 

O(E) 

+ [tree level contribution from the NLO Lagrangian] 
s >, ' 

0(£2) 

+ [one-loop contribution from the LO Lagrangian] 

^ v ' 

0(£ 3 ) 

+ [ tree level contribution from the NNLO Lagrangian] + 0(E 4 ), (26) 

0(E 3 ) 

where [i c is the vector meson mass in the chiral limit. In fact, the second term on the right hand 
side of eq. 1)26(1 comes from the 0(a) Lagrangian only. As has been shown in Ref. the 0(p) 
term does not contribute at tree level and the chiral correction to the vector meson mass starts at 
0(E 2 ). 

In this paper we are only interested in deriving the vector meson masses through 0(E 3 ). In that 
case we neither need the NLO terms of 0(ap,p 2 ) nor the NNLO terms of 0(mp, p 3 , ap 2 , a 2 p) . The 
reason is that all these terms contain at least one momentum factor p. By using the lowest order 
equation of motion, we can replace the residual vector meson momentum with the momentum p„ 
of the pseudo scalar and a contribution to proportional to a. 3 The resulting term has a part with 
at least one p n and a part of 0(a 2 , am, a 3 ). For the former, p w ~ d^ir implies that this term always 
contains at least one pseudo scalar field together with a derivative d^. This kind of term results in a 
three-point vertex (vector-vector-pseudo scalar) as a non-vanishing leading term and does not give 
a tree level contribution to the vector meson mass. Hence the terms of 0(ap,p 2 ,mp,p 3 ,ap 2 ,a 2 p) 
are not needed for our calculation and we do not list these terms in the next subsection. 



C. Effective Lagrangian 

In this section, we show our results for the terms in the chiral Lagrangian which contribute 
to the vector meson mass through order E 3 . By construction it involves the meson fields and is 
invariant under Lorentz and chiral transformations, charge conjugation and parity. 4 

0(p): 



3 The leading order equation of motion, (—iv ■ d + aa s )5 M , is easily read off from the LO chiral Lagrangian, see eq. 
\A2l in appendix I A II 

4 Just for convenience we present our results in Minkowski space, which does not make a difference for the com- 
putation of the vector meson masses. The underlying lattice theory, however, is usually formulated for Euclidean 
space time. 
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C p = -iStivd)Si* 

(27) 



-i{Ol{v-V)0") 



+ig x (sl(O v A x ) - S^OlAxYjv^** 
+ig 2 ({OlO u }A x )v a e^, 

with (X) = tr(X) for flavor indices. The derivative id^ acting on the vector mesons provides 
the residual momentum r^ = k fl — fiv v u- We therefore consider the kinetic terms for the vector 
mesons as 0(p)- The terms proportional to g\ and g 2 involve A\, defined in eq. ()23j) . Expanding 
in terms of the pseudo scalar fields this quantity starts as A\ = <9\il//+0(n 3 ). The presence of 
the derivative means that both the g\ and g 2 term are of 0(p). 

As is well known, the vector mesons are not stable. This effect can be taken into account by 
including an anti-hermitian term in the Lagrangian whose coefficient is proportional to the decay 
width. However, the decay width is rather small and this contribution is usually ignored We 
also neglect this term in this paper. 

0(a): 

L a = ai (W + )SlS» 

+a 2 ((Olw + )S» + (0,W + )S^) 

(28) 

+a 3 ({OlO»}W + ) 
+au{W+){OlOP). 

This 0(a) Lagrangian is a new element of this study. Since the spurion field A associated with the 
lattice spacing transforms exactly like the quark mass spurion field (cf. eq. p9[l). this part of the 
Lagrangian has the same structure as the following 0(m) Lagrangian. 
0(m): 

C m = X s (M + )slS" 

+\ os {(Oj 1 M+)S» + {O^M + )S^) 

+\ ol (M + )(OlO») (29) 

+\o2({OUO»}M + ) 

+A f j,sLs». 

Note that the leading order mass terms /j, s Sj l S fl and ^(Oj^O^), where pt s and // Q are the masses 
for the singlet and octet vector meson in the chiral limit, respectively, are absent in the 0(p) 
Lagrangian. Their effect is already included through the phase factor exp(— irrtyv ■ x) in the 
definition of the heavy meson fields. The mass difference 

A/i = [j, s - ii Q (30) 
fO 



between the singlet and octet vector meson mass, however, needs to be introduced explicitly. This 
term, which is the last term in eq. Q29j). is considered to be of O(m), since A/i < 200 MeV ~ m s 

0(a 2 ) : 

C a2 = faiW+W+^lS" 
+ (3 2 (W + ) 2 sis» 

+ 3 {{O} l W+W+)S'> + {O^W + W + )Srt) 
+ fo(W + )((0l l W+)S>* + {O^W + )S^) 

+p 5 ({ol,o»}w + w + ) 

+ (3 6 (Oiw + 0»W + ) (31) 
+ [3 7 {W + )({Ot,0»}W + ) 

+ P 9 {W + ) 2 {Ol<D») 
+ (3 w (Oiw + )(0»W + ) 
+ {W+ W-). 

In total we find twenty terms, however, the parts involving W_ are not needed. By expanding in 
powers of the pseudo scalar field one finds W- = (2ia/ f) H at leading order. Hence, all these terms 
have at least two pseudo scalar fields and therefore do not contribute to the vector meson masses 
at tree level. 

5 We assumed the definition my — n a . If we had chosen my = Us the last term in eq. (1291 would be — A/^CJtC*). 
Our final results for the vector meson masses are, of course, independent of this particular choice. 
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0{am): 



Cam = ll{W + M + )S\ L S^ 

+ l2 (W+)(M + )SLS» 

+ 73 (M + )((Olw + )S» + {O^W + )S^) 

+ 74(4<0"{M +) W + }) +S tl (Ort{M + ,W+})) 

+ l5 (W+)((OlM + )S» + (O^M + )S^) 

+ l6 (M + )({Olo^W + ) 

+ yr(OlO' i )(M + W+) (32) 

+ l9 ({Ol,0^{M + ,W + }) 

+ 7i «Oj,M+O'*W + ) + 

+ 7ii(^ + )({0i,On^+) 

+ 7i 2 ((0) t M+)(^W / +) + (0^+)(0 M M+)) 

+ (M + ,W+ -> M_,W_). 



One can easily check that M_ ~ II at leading order. Therefore, the terms involving M_ and W- 
do not contribute to the vector meson mass at the order we are working. 

For our purposes it is not necessary to list the full O (a 3 ) La grangian, which is quite cumbersome. 
Since we will not encounter any divergences proportional to a 3 in our calculation, we do not need 
the coefficients in the 0(a 3 ) Lagrangian as counterterms. The tree level contribution of the 0(a 3 ) 
Lagrangian merely gives an analytic a 3 term to the singlet and octet vector meson mass, which 
we can simply add to our final one loop result. Note, however, that the 0(a 3 ) Lagrangian does 
not give an off diagonal contribution to the singlet-octet two point function at 0(E 3 ). All possible 
ways to take the trace in terms like 



yield at least two pseudo scalar fields in the non-vanishing leading term when one expands W± in 
powers of the II field. By noting (O^) = 0, W+ = 1 + 0(n) and W- = 0(ir), one can easily check 
this statement. 

Having determined the chiral effective Lagrangian we can set the spurion fields to their constant 
values, i.e. we take A — ► al in eqn. ((2"T|) - (|2"§)) and eqn. (|3*Tj) - (f3*2"|) . This replacement simplifies 
the fields W± to 



S\ Jl O tl W+W + W+ + /i.e. or 



SjCW + W_W_ + h.c. 



(33) 



W ± - |(E±Et). 



(34) 
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The Lagrangian in terms of the component fields 7r a , p® for the pseudo scalar and vector mesons 
is presented in appendix [XJ together with the Feynman rules relevant for the calculation of the 
vector meson masses. 



D. Vector meson masses 



The full propagator for the p and the K* is given by 

iA ab (k) = ~ %PiiV 6 ab + ~ iP(ia [-iP^Sffffc)] ~ lPpu +••• 
h v ■ r — aa Q i v ■ r — aa Q \ v ■ r — aa \ 

= -iPfiu[v ■ k - p ~ aa Q i + Eo(fc)]^ 1 , (35) 

where a, b runs from 1 to 7. Recall that the external four-momentum for the octet vector meson is 
written as k v = p Q v u + r v where p Q is the octet mass in the chiral limit. Using the Feynman rules 
summarized in appendix 1X1 we obtain for the octet-to-octet self energy through 0(E 3 ) 

So (*) = -[(NMoKi + a 2 (3 Q + NaM Qlol )5 ab + (A o2 + a lo2 )(2M 6 aa + M 8 d ab8 )] 



+^5 ab 2C 21 (w s ,m a ) + M Y,d ace d bce 2C 21 (w Q ,m c ) 

J c,e 

+ a Tf2 Sab E J o(-c) + ^ E d abe d-I (m c ), (36) 



where 



w ,s = v ■ r - aa ol)S , (37) 

and d abc is the totally symmetric d-symbol of su(iV). 6 The coefficients a Q i are particular combi- 
nations of the low energy constants a\, . . . ,«4 which enter eq. (|28|). Similarly, the coefficients (5 Q 
and 7oj are combinations of the low energy constants Pi, ... , /3io and 71, ... , 712, respectively. The 
full expressions for these combinations can be found in appendix even though it will not be 
necessary to keep track of the original constants in the following. The functions C 2 i(lo,M) and 



6 We keep the number of flavors N undetermined at intermediate stages since N provides some useful checks for our 
equations. In the final results we will set N equal to 3. Note, however, that intermediate results like eq. (|3fi^ do 
not hold for arbitrary N unless all masses are degenerate. 
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Io(M) stem from the loop integral and are denned as 

C 21 (w,M) = \[{M 2 -w 2 )J(w,M)+wI (M)] 



w 



'M 2 

127T 2 I 2 



w 
~3 



w 



J(w,M) = 



R + \n- 



M 2 



+ 



1 



I (M) 
R 



8tt 2 

M 2 
16vr 2 
2 



2 \/ M 2 — ic 2 arccos 



M 



for w 2 < M 2 , 



[ln(4 7 r)+r / (l) + l], 



(38) 

(39) 

(40) 
(41) 



n — 4 

with the dimension n. The one-loop integrals encountered here are essentially the same as those 
in baryon chir al p erturbation theory, and a detailed derivation of the formulae (|38j) - (|41j) can be 
found in Ref. [12]. 

In order to cancel the divergence in R, defined in eq. (|41|) . the O(am) coefficients 7j (i = ol, o2) 
need to be properly renormalized, 



7t = 7» (A 1 ) + —^kR, 



(42) 



8^ 2 / 2 ' 

where the coefficient B is a low-energy parameter in the LO chiral Lagrangian for the pseudo 
scalars (see Ref. []] and appendixEJ). The coefficients 7[(/x) denote renormalized parameters, fi is 
the renormalization scale 7 , and the coefficients Cj are given by 



Col 
Co2 



iV 2 



1 



a Q i- 



iV 



«o2 



3iV : 



6iV' 



(43) 
(44) 



The on-shell condition for the p and K* meson (a = 1, • • • , 7) is given by (no sum over a) 



■k- Ho- aa Q i + So (&), 



ki. 



m a v 



k 1 



(45) 



where m a is the physical mass. So far we have chosen to express the self energy as a function of the 
four-momentum ku- ^o'{k) may also be parametrized by the scalar variables w Q and w s , defined in 
eq. (|37j). and we can write S a (A;) = S a (u; , w s ). It is easy to see that the definition of w Q together 
with the on-shell condition (|45|) yields, 



Wr 



m a - Ho - aa Q i 



-^{k = m a v). 



(46) 



In the following we will simply write ^ for the renormalized parameters and suppress the dependence on /x. 
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w is considered to be of 0(E ' ) since the self energy Eq 1 is of 0(E). For w s we find the similar 
result 

w s = m a - Ho - aa s = -E °(A; = m a v) - a(a s - a Q i). (47) 

Here we assume that the difference between the leading lattice artifacts for the octet and the singlet 
mass is of order E 2 and therefore small 8 , 

a(a s - Q i) = aAa = 0(E 2 ). (48) 

This is supported by a strong coupling analysis performed in Ref. [l^]. Expanding the effective 
potential through second order in terms of the vector meson fields, it has been shown that there 
exists no difference between the singlet and octet vector meson to all orders in the strong coupling 
expansion. Whether w s and w Q is small away from the strong coupling limit can and should be 
checked in the actual lattice simulation. 

Provided the assumption (|48|). both w s and w are of 0(E 2 ) and we can set io 0)S = in E Q a , 
since the difference is beyond the order we are working to. We therefore conclude that 

m a = /u + aa Ql - E °(«; = 0, w B = 0). (49) 

Note that the function C21 becomes much simpler when the first argument is set to zero, 

C 21 (w = 0,M) = ^-M\ (50) 

247T 

Let us now consider the sector of the p s and the (jp mesons. Due to mixing, the propagator for 
this sector is written as a 2x2 matrix, 

1 



v ■ r - aa ol + T, (k) S s(fe) 

Eos (k) v ■ r - aa s + E s (k) 



(51) 



where the self energies for octet-to-octet (a = 8), singlet-to-singlet and octet-to-singlet are given 



If this assumption is not valid, then w s is of O(a) ~ O(E). In this case, the function C2i[w a ,rha) in eg. 13611 
contains divergences whose coefficients are of order am and a 3 . Such divergences would alter the renormalization 
conditions for 7 i, 2 in eq. 142H . Also the 0{a i ) counterterms may receive a divergent renormalization in order to 
cancel additional divergences. Furthermore, if w B is of 0(a), the condition w 2 < M 2 in eq. (1391 may not be valid. 
Hence the functional form of J(w, M) might be changed according to the relation between w and M as in Ref. 
eq.(C27). 
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by 



X a = 8 > b=8 (k), 



E s (fc) = -(NM \ s + a 2 (3 s + NaM ~f s + A/x) 



2a K 



(52) 
(53) 



£os(*0 



(A os + a7os) 



+ 



.92 



^d cc8 2C 2 i(w ,m c ) + 



a, 



3 >/2 



^d cc8 / (m c 



(54) 



For the cancellation of R, in a similar way to the case of p-K* sector, O(am) coefficients ji (i = s, os 
should be renormalized as in ea. ()42j) with 

N 2 -l 



C B 

Cos 



N 2 ' 

5 

; 6iV' 



(55) 
(56) 



dct 



0. 



(57) 



The physical masses mi, i = 1,2, for this sector are determined by the on-shell condition 

v ■ r - aa Q i + T, (k) E s(&) 

Sos(^) v ■ r - aa s + E s (fc) 

at kfj, = rriiV^. The task is to solve the equation 

m — /jl — aa ± + Eo(fc = mv) Eos(^ = mv ) 

Sos(^ = mv ) m ~ Ms ~~ aa s + E g (/c = mv) 

with respect to the mass rrii, where we have defined E g as Eg = E g — A/i. The solutions of the 
equation are given by 

1 



det 



0. 



(58) 



m± 



fi s + fi Q + a(a s + a ol ) - (Eg + S Q ) 



±y/{A» + aAa - (E' s - Sq)}" + 4(E OS ) 2 



(59) 



Assuming that the SU(3) breaking (m — m s )/3 is smaller than the average (2m + m s )/3 of the 
quark masses 9 together with {Afj, + aAa — (Eg — Eo)} > 0, the physical masses are given by 

(Eos) 2 



m+ = /i s + aa s - Eg + 
= Ho + aa \ - E ■ 



m 



A/i + aAa - (E' g - E ) ' 

(Eos) 2 
A/i + aAa - (E' g - E Q ) 



(60) 
(61) 



9 This assumption means that the diagonal part is dominant compared to the off-diagonal element, {Afi + aAa 
(E^-Eo)} 2 >4(E OS ) 2 



16 



As in the case of the p and K* sector, we consider the self energies So, OS an d S g as functions oiw 0)S 
when we impose the on-shell condition. Setting the physical masses equal to m±, the parameters 
w os are given as 



w = m — fi Q — aa \ 

(Sqs) 2 



Afj, + aAa - % - AM+aAa ° s ( ^_ So) > for m = ?n+, 
~ S o- A M +a£-(sVso) ' for m = m_, 



w s = m — p Q — aa ; 



s 



A/x-£' g - A/ ,+aA«- S (4-£o) ' for m = m + , 



^2 



-aAa - Sq - A^l^fw So) , for m = m_. 



(62) 



(63) 



Hence it turns out that both w Q and w s are of 0(E 2 ), since A//, aAa, So,s an d Sos ar e of 0(E 2 ). 
Therefore, we can set w QiS = in the self energies, as we did in the p-K* sector. 

With these preparations our one-loop results for the vector meson masses are as follows: 



17 



m. 



m + 



m_ 



m 



(88) 



™(00) 



m 



(08) 



m Q (a) + X x (a)x + 2X y (a)y 



1 



(5? + gfli)(a: + 2y)5 + 2g\{x - y)i + -&{x - 2y)5 



12tt/ 2 
2a 

mo (a) + A z (a)a; - A y (a)y 



3a Q i + a o2 ) L,,- + ( 4a D i - -a o2 ) La + ( a i - -a 2 ) L 



1 



12tt/ 2 
2a 
3p 



2' 
3a Q i 



+ 2y)2 + 5 ? + - 5 1 (x - y)2 + -g 2 2 ( x _ 2 y)2 



1 



1 



«o2 ) L n + 4a Q i + -a Q 2 £ftr + a i + -a Q 2 £ 



1 



6 



■/?? 



(08) 



m (00) ~ m (88) 
2 



^(00) + 

TO(88) - 

m (00) ~ m (88) 

mo (a) + X x (a)x - 2\ y (a)y 



(08) 



1 



2g\{x + 2y)i+ -g\(x - y)? + ( g\ + -g\ [x - 2y) 



12tt/ 2 
2a 
3p 

m s (a) + a x (a)x 



3a i - a Q 2 ) L n + ( 4a D i + -a o2 ) La + ( a i + ^a 2 ) L 



1 



12vr/ 
2a 
3/2 

(7j,(a)y 

12^P 9192 U 
a 



3 3 3 

3(x + 2y)2 + 4(x - y)5 + (x - 2y) 2 



3L ff + 4La + I,, 



2r 



3 3 3 

3(x + 2y)2 - 2{x - y)2 - (x - 2y)z 



3/ s 



3L n — 2La — 



2m 



(08) 



tan 20 = 

m (oo) ~ m (88) ' 

where 6 denotes the mixing angle in the p 8 -</>° sector. We introduced the two parameters 

2B 



x 



-(2m + rh s ), 



D 



y = —{m — m s 



as a short hand notation for convenient combinations of quark masses [l| . The leading order pseudo 



Q. 



(64) 



(65) 
(66) 
(67) 



(68) 



(69) 



(70) 
(71) 



(72) 
(73) 
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(74) 




scalar masses assume a very simple form in terms of x and y, 

fh% = x + 2y, 

rh 2 K = x-y, 
m 2 = x — 2y . 

For the chiral logarithms stemming from the loop integration we have introduced 

~ 2 / ~ 2 ' 

^ = SM^-). ( 75 ) 

The functions 

mo(o) = + 14a + F 2 a 2 + T/ 3 a 3 , (76) 

ms(o) = S'o + Sia + ^a 2 + 5 , 3 a 3 , (77) 

A*(a) = Af + A«a, A. y (a) = Af + A«a, (78) 

o-«(a) = of + a^a, o- tf (o) = of + o^a, (79) 

parametrize the analytic lattice spacing dependence. We have, for convenience and for the sake 

(i) (i) 

of transparency, introduced new constants Vi,Si, i = 1, . . .4 and \x,y, &x,y, j = 1,2. These new 
constants are combinations of the low energy constants in the chiral Lagrangian. They are explicitly 
given 

(80) 
(81) 
(82) 

(83) 

(84) 

(85) 



V 


= fJ'o, Vi = a \, V 2 


= A>, 


So 


= Ms, Si = a s, 5*2 = 


= &, 


Af 


N\ i + 2A 2 . (i) 
2B 


^7o r i + 27o r 2 
25 


Af 


Ao2 ,(1) 7o2 

B ' y B ' 






^A s (1) iV 7s r 
25 ' x 2£ ' 




4 0) 


B ' CT f B 





However, from a practical point of view there is no need to keep track of the original low energy 
constants. Note that we have added an 0(a 3 ) term for the octet and singlet vector meson mass, 
but not for the singlet-octet mixed contribution ?7i(08), as discussed in sect. IIII 
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Note that there are chiral logarithms present in the results for the vector meson masses. These 
chiral corrections are a lattice artifact, as can be seen from the presence of the a in front of the 
logarithms. This chiral log contribution vanishes in the continuum limit and the mass formulae 
converge to the results given in Ref. However, as long as the lattice spacing is non-zero, there 
are two kinds of non-analytical quark mass dependence, proportional to m^ 2 and m q log m q . It will 
be very interesting to study the competition of these two contributions in actual lattice simulations. 

Let us count the number of unknown parameters. In the case of the p and K* meson, we find 

(i) (i) 

twelve unknown combinations of low-energy constants: Vi, Ai , Ay , g\, g2, ol \, «o2- However, 

as long as we are interested in performing fits to lattice data at a given and fixed lattice spacing 

a, this number essentially reduces to seven independent constants: too, X x , \, <7i> 92i a oi> Oo2- 

(i) (i) 

For the m± sector, we have additional ten low energy constants: Si, Ox, a\ , , a os , a s . When we 
consider at a fixed lattice spacing, there are essentially five independent constants: ms, cr x , a y , 



IV. CONCLUDING REMARKS 



We have derived the one loop expressions for vector meson masses using an effective theory 
based on the heavy vector meson formalism of Ref. Q] . The effects due to a non- vanishing lattice 
spacing introduce a fair number of new unknown low-energy constants in the chiral Lagrangian. 
However, the actual number of fit parameters in the mass formulae seems small enough for the 
expressions to be useful for the chiral extrapolation of actual lattice data. In the case of the p and 
K* meson mass, the number of unknown fit parameters is seven. The numerical simulations of the 
CP-PACS/JLQCD collaboration are carried out with five different (degenerate) up and down quark 
masses and two different values for the strange quark mass. The total number of independent data 
points is 18 without data for the <j> and uj mesons. If one measures these masses including noisy 
disconnected contributions, the number of data points increases to 36. In both cases the number 
of data points is well above 7. 

Here we have calculated the expressions for the vector meson masses only. In the calculation 
of the vector meson decay constants some additional low-energy constants enter, depending on 
the choice for the vector current. Apart from this the calculation is straightforward and currently 
under way j^ . 

Our main motivation for the calculation presented here was to derive fit forms for the chiral 
extrapolation of the vector meson masses. If our formulae describe the lattice QCD data well, 
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we obtain, as a byproduct, estimates for the low energy constants of the chiral effective theory for 
vector mesons. Particularly promising are the leading order parameters g\ and gi- These couplings 
are not corrected by lattice artifacts through 0(E 3 ). Provided that the corrections of 0(E 4 ) and 
higher are negligible, we have a good chance to determine g\ a nd a ? in a fit, even if lattice data is 
available for one lattice spacing only. The chiral quark model [2l| predicts g 2 = 0.75 in the large 
N c limit where g\ = -^92- On the other hand, the estimate g 2 ~ 0.6 is given in Ref. Q, where 
vector meson XPT is applied to r decay processes together with a comparison of the results with 
experimental data. If we are able to extract g\ and gi in a fit to lattice QCD data, we obtain 
results based on first principles QCD without relying on model dependent assumptions or the large 
N c limit. 

The heavy vector meson formalism of Ref. P] is not the only way to derive an effective theory 
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for vector mesons. Another approach employs the so-called hidden local symmetry 
approach has a definite counting scheme, and in principle one can formulate it for non-zero lattice 
spacing too. It would be interesting to perform the calculation of the vector meson masses based 
on this alternative approach. The comments in the previous paragraph about the extraction of low 
energy constants through a fit to lattice data applies to this effective theory as well. 
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APPENDIX A: FEYNMAN RULES 



In this appendix, we present the Feynman rules needed in the calculation for the vector meson 
masses. The lattice spacing a shown in the following equations stems from the spurion field W± 
once the replacement A — ► al has been made. 
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1. Expanding the Lagrangian 



In terms of the component fields, the LO Lagrangian which includes up to two pseudo scalar 
meson fields is expressed as 

= (£g + O + C 1 / + (Cf + CI-) + 0(tt 3 ), (Al) 
C^ + C^ = Sl{-ivd + aa s )S» + Y J P a J{-™-d + aa ol )p a » 1 (A2) 

a 

C l n = ^E(4^-^^ t )^^e^ ACT + ^1^46^V^A7r c ^e^, (A3) 

•* a J a,b,c 

Cf + Cf = -4 E fadehcep a ^ b (vd* C )p d » 
a,b,c,d,e 

a,b,c 
2 

a,b,c,d e 

- a Wf2 E SabSlSWv", (A4) 
where the newly introduced coefficients relate to the ones defined in eq. (|28[) according to 

a s = iVai, (A5) 

a os = Not2, (A6) 

Ooi = 2a 3 + iVa4, (A7) 

a o2 = Ar a3 . (A8) 

The summation over roman indices is understood to range from 1 to 8. f a b c and d a b c are the 
structure constants and the totally symmetric d-symbol of su(N) respectively. Obviously, we are 
most interested in the case N = 3. 

The NLO and NNLO Lagrangian results in the following tree level contribution to the vector 
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meson masses: 



£NLO + ^NNLO = (£m + £ a 2 ) + £am 

= C^ + C^ + C^ + 0(ir), (A9) 
£m+£°J + £°am = (N M \ s + a 2 (3 S + N M a ls + A^S^ 

+ (A os + a 7 os)^ £ f«8(#S" + 

+(iVAf Aoi + a 2 /3 D + iVM «7oi) ^ p£V a 

a 

+(A o2 + a7o 2 ) ^(2M <5a6 + Mfid <M )p$P>* i (A10) 



where 



A = Nfa + N 2 ^, (All) 

fl, = 2/3 5 + /3 6 + 2Ar/3 7 + Ar/3 8 + Ar 2 /3 9 , (A12) 

7s = 7i + A r 72, (A13) 

7os = 274 + ^75, (A14) 

701 = 27 6 + 7 7 + A^7 8 , (A15) 

702 = 279 + 710 + A r 7n- (A16) 



2. Propagators 



The propagator for the singlet and octet vector meson is given by 

-iP u 



GUr) = ^ , (A17) 

^ v ■ r — aa s 

G°f(r) = ~ iP ^ , (A18) 
r v ■ r — aa Q i 



where is the residual momentum. 



For the pseudo scalar meson, the propagator is 



G?S(P) = "1^2 > ( Al9 ) 
— m„ 



where the leading order pseudo scalar mass fh a of Ref. [1|| is used. The explicit form of the mass 
is 

m 2 a = m 2 (l - Nc 3 - c 3 ) - m 2 v 7Vc 3 - Nc 2 - c 2 , (A20) 
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with 



m: 



N2 _ l ^2 m2 a> 



(A21) 



where m\ is the meson mass in the continuum limit (i.e. m\ = m| = m\ = mf\, vnh = m\ = m\ 



m 7 = m ?o m 8 = m ri)- c i' c 3> ^3 are 0( a ) terms Ql, and C2, C2 are 0(a 2 ) terms 



y. 



They are 



low energy constants in the pure pseudo scalar meson sector. The mass m\ is expressed in terms 



of the shifted quark masses, 



m n 



ml = 2Bm, a = 1,2,3, 

m\ = B(m + m s ), a = 4, 5,6,7, 
mi 



v 



(A22) 



9 25 . , 
- -^-(^ + 2m s ), a 



Vertices 



NLO vector-vector 2-point vertex from £j^ r + + £ 
v, b /Li, a 



-(2M S ab + M 8 d ab8 )(A o2 + a7 o2 )], 



> X > iP^(NM X 3 + a 2 (3 s + NaM j s + A/x). 



The oriented narrow line represents the singlet vector meson and the oriented wide line 
represents the octet vector meson. 
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LO vector-vector-pseudo-scalar 3-point vertex from 

p, b 



v, a 



l -fp x v a e^ x ° 5 ab , 



P, b 



Mj a 



ifp x v a e^5 ab 



P, c 



v, b 



H, a 



^pxv„e^ x °d abc 



The dotted line represents the incoming pseudo scalar meson. 



LO vector-vector-pseudo scalar-pseudo scalar 4-point vertex from C 



2tt 

a 



d 



v, b 



\ * 
\ t 
\ / 
\ t 
s t 
s t 



v, a 



\ * 
\ t 
\ / 
\ t 
s t 
s t 



\ I 
\ t 
S t 
\ t 
\ t 
\ / 

> Y > 



25 



[1] S. Aoki, 0. Bar, T. Ishikawa and S. Takeda, |hep-la t/0509049 

[2] S. Aoki, O. Bar, T. Ishikawa and S. Takeda, work in progress. 

[3] T. Ishikawa et al, Nucl. Phys. Proc. Suppl. 140 (2005) 225. 

[4] S. R. Sharpe and J. Singleton, Robert, Phys. Rev. D58 (1998) 074501. 

[5] W.-J. Lee and S. R. Sharpe, Phys. Rev. D60 (1999) 114503. 

[6] O. Bar, Nucl. Phys. Proc. Suppl. 140 (2005) 106. 

[7] H. R. Grigoryan and A. W. Thomas, |hep-lat/0507028 



[8] E. Jenkins, A. V. Manohar and M. B. Wise, Phys. Rev. Lett. 75 (1995) 2272. 

[9] S. Aoki, Phys. Rev. D68 (2003) 054508. 

[10] Y. Namekawa et al., Phys. Rev. D70 (2004) 074503. 

[11] E. Jenkins and A. V. Manohar, Phys. Lett. B255 (1991) 558. 

[12] J. Bijnens and P. Gosdzinsky, Phys. Lett. B388 (1996) 203. 

[13] J. Bijnens, P. Gosdzinsky and P. Talavera, Phys. Lett. B429 (1998) 111. 

[14] J. Bijnens, P. Gosdzinsky and P. Talavera, Nucl. Phys. B501 (1997) 495. 

[15] M. Booth, G. Chiladzc and A. F. Falk, Phys. Rev. D55 (1997) 3092. 

[16] C.-K. Chow and S.-J. Rey, Nucl. Phys. B528 (1998) 303. 

[17] G. Rupak and N. Shoresh, Phys. Rev. D66 (2002) 054503. 

[18] S. Scherer, Advances in Nuclear Physics. Vol. 27 (2002) 277. 

[19] S. Aoki, Phys. Rev. D34 (1986) 3170. 

[20] S. Aoki and S. Takeda, work in progress. 

[21] A. Manohar and H. Gcorgi, Nucl. Phys. B234 (1984) 189. 

[22] H. Davoudiasl and M. B. Wise, Phys. Rev. D53 (1996) 2523. 

[23] M. Harada and K. Yamawaki, Phys. Rept. 381 (2003) 1. 

[24] O. Bar, G. Rupak and N. Shoresh, Phys. Rev. D70 (2004) 034508. 



26 



